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Abstract 

We discuss the dynamical analysis in f{R,T) gravity (where R is 
Ricci scalar and T is trace of energy momentum tensor) for gravitat¬ 
ing sources carrying axial symmetry. The self gravitating system is 
taken to be anisotropic and line element describes axially symmet¬ 
ric geometry avoiding rotation about symmetry axis and meridional 
motions (zero vorticity case). The modified field equations for axial 
symmetry in f{R, T) theory are formulated, together with the dy¬ 
namical equations. Linearly perturbed dynamical equations lead to 
the evolution equation carrying adiabatic index T that defines impact 
of non-minimal matter to geometry coupling on range of instability 
for Newtonian (N) and post-Newtonian (pN) approximations. 

Keywords: f{R,T) gravity; Axial symmetry; Instability range; Adia¬ 
batic index. 


1 Introduction 

Recent developments in astrophysics and structure formation theories reveals 
that gravitating sources might deviate from most commonly studied spheri¬ 
cal symmetry. Such deviations in realistic scenarios appear incidently, giving 
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rise to the worth of non-spherical symmetries in gravitating objects. Herein, 
we intend to look into the implications of restricted class of axially symmetric 
sources (avoiding reflection and rotation) on gravitational evolution in con¬ 
text of f{R, T) theory of gravity. Consideration of dynamic sources together 
with the angular momentum is a cumbersome task, however, observational 
data suggests that the lack of spherical symmetry prevails the more practical 
and worthwhile situations. A viable f{R,T) model (^ > 0, ^ > 0) with 
locally anisotropic matter distribution has been taken into account for the 
dynamical analysis. 

The evolution of gravitating sources has been studied with a great deal 
of interest in recent past. Stars tend to collapse when outward drawn pres¬ 
sure decreases because of continuous fuel consumption, leading to imbalance 
in outward forces and inward acting gravitational pull. In such situation, 
gravitational force becomes the only governing force, massive stars burn nu¬ 
clear fuel more rapidly and so more unstable as compare to the stars with 
relatively less mass. There are many other factors other than mass of the 
gravitating source that implicate intense modihcations in range of stabil- 
ity/instability such as isotropy, anisotropy, shear, dissipation and radiation. 
Chanderashekar [1] shared valuable explorations to set instability range for 
spherically symmetric gravitating source in the form of adiabatic index T 
comprising pressure to density ratio with time transition. 

Hillebrandt and Steinmetz [2] presented the instability criterion for anisotropic 
matter conhguration of gravitating objects. Herrera et ah 0-111 contributed 
majorally in establishment of instability range of general relativistic fluids 
for different cases (isotropic, anisotropic, dissipative collapse etc), remarked 
that pressure anisotropy largely participates in setting dynamical instability. 
Moreover, they also worked out the imprints of axially and reflection sym¬ 
metric static and dynamic sources by a general framework and some analytic 
models [8]. Axially symmetric shearing geodesic and shearfree dissipative 
fluids are discussed in [9l [10], where shearing geodesic case represents the 
zero radiation production. 

General Relativity (GR) is a self-consistent theory, it is adequate for the 
explanation of many gravitational phenomenons up to cosmological scales. 

The scheme of GR appears to disagree with progressing observational data 
such as large scale structures ranging from galaxies to galaxy clustering, 

I A-type Supernovae, cosmic microwave background nu-mi etc. Alterna¬ 
tively, it can be said that GR is not the only dehnite gravitational theory 
that is suitable for all scales. Many attempts have been made to validate 
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gravitational theories on large scales and coup with the cosmic acceleration 
[IS]-[21], by introducing modified theories of gravity [22]-[30] for e.g. f{R), 
f{G), Brans-Dicke theory, f{R,T) and so on. 

Since the introduction of f{R,T) theory in 2011 [31], people [32]-[M] 
worked on energy conditions along with its cosmological and thermodynamic 
implications. The f{R,T) theory represents generalization of f{R) theory 
carrying non-minimal matter to geometry coupling. Extensive work has been 
done on instability range of spherically symmetric stars in GR as well as in 
modified theories of gravity. Literature on dynamical analysis of axially 
symmetric sources can be witnessed in GR. However, being a heavier task to 
handle modified dynamical equations in modified theories, very few attempts 
have been made to explore axial symmetry. 

The purpose of this manuscript is to work out the instability problem for 
axially symmetric (in absence of reflection and radiation) anisotropic sources 
in context of f{R,T) gravity. The reason of avoiding reflection and rotation 
terms in axially symmetry is only to somehow reduce the complications in 
analysis. The modified EH action in f{R,T) admits arbitrary function of R 
and T to account the exotic matter. The action in f{R,T) is given by [31] 

+ ( 1 - 1 ) 

where £(m) represent matter Lagrangian and g represents the metric. Several 
choices of £(m) can be considered, each of which stands for a specific form of 
fluid. 

The article arrangement is: The matter configuration and components of 
field equations together with the dynamical equations are furnished in section 
2 . Section 3 covers the information about the f{R,T) model and perturbed 
conservation equations leading to the collapse equation. Section 4 contains 
range of stability of N and pN limits in the form of adiabatic index. The last 
section consists of concluding remarks followed by an appendix. 

2 Interior Spacetime and Dynamical Equa¬ 
tions 

The general line element for axially symmetric compact objects constituting 
five independent metric coefficients is given by 

-|- B'^dr^ + B^r^dO^ + C'^dcp'^ + 2Gdtd6 + 2Hdtd(j)^ (2.2) 
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where the metric functions A, B, C, G, H have dependence on time, radial and 
axial coordinates {t,r,6). Herein, we have ignored the meridional motions 
and rotation about the symmetry axis. Absence of dtdd and dtdcf) terms lead 
to the restricted character i.e. vorticity free case. The modihed equations 
are highly non-linear in nature, so it is a tough task to handle such equations 
with non diagonal entries in the metric tensor, that is why we have taken zero 
vorticity case to somehow manage dynamical analysis by analytic approach. 

The reduced form of general axial symmetry with three independent met¬ 
ric functions is [35] 

ds^ = r, 6)dt^ + r, 6){dr‘^ + r^dO"^) -|- r, 6)d(j)^. (2.3) 

Taking £(m) = —p, SttG = 1 and varying the EH action fl3.28p with respect to 
the metric tensor puv leads to the following form for modihed held equations 

(/t + + PQuvfr H ^- Quv 

+ (VuVu — guvG\)fB\ 1 (2-4) 

where □ = V“V„, fn = df{R,T)/dR, J't = df{R,T)/dT, is covariant 
derivative and is the energy momentum tensor for usual matter. The 
matter conhguration is considered to be locally anisotropic [2H], given by 

= {p + P±)VuV, - {K^K, - - P,,) - (L,L, - 

— Pxx)-\-Pguv + ‘^R{uLv)Pxy, (2-5) 



where p is the energy density and 


P 


+ Pyy + Pzz)') 


Quv KK, 


Pxxi Pyy, Pzz ^ud P^y are respective stresses causing pressure anisotropy, pro¬ 
vided that Pxy = Pyx and Pxx 7^ Pyy 7^ Pzz- Ku and Lu represents the four 
vectors in radial and axial directions, respectively and Vu is for four-velocity, 
these quantities are linked as 


Vu = -A6l K^ = B5lG Lu = rB5l. ( 2 . 6 ) 
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The components of modified (effective) Einstein tensor are 


^00 


G02 

Gii 


G^^ 

G22 


G33 


A^fR 


p + 


fe 

Jr 


A^fR 

2B^ 


f-RfR Jr h'B ^ C\ fj, fl , 2B' 


— I-h— 1- — I- + 

A^\BG' 1 


52 \r 


B 


J.2 J^2 
-1 

-1 

r^A^B'^fn 
1 


C^\,fR 

G B^ 


B 

~^fR + (g/fi “ fn 


^fR + ^Ir - fR 


B^fR 


r. 


'^2 ir A^ B G 


Pxxifr + 1 ) + pfr + 


A' B' 


Pxyifr + 1) + /ij 


r^B^fR [ 
1 

r^B^fR 



f-PfR /. 


Gd 

~C 


^’„(/t+1)+/>/t + ^ +^-- 


C" 


(2.7) 

( 2 . 8 ) 

(2.9) 
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( 2 . 10 ) 
( 2 . 11 ) 
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(2.13) 


Here dot, prime and 6 indicate the derivatives w.r.t t,r and 9 coordinates 
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respectively. The expression for Ricci scalar is 


R = 


2 A 12B CA B B 2C\ 2B C 

R 2 a[^^cI~b[b^~cI~^~c 


2 \A" A'C B" \ (A! B' C 


B"^ [A AC B r V ^ B C J B^ 
C 1 fA^^ B^^ R® 2 A^C^y 


1 fA^^ B^^ 


(2.14) 


In order to explore stellar evolution, one needs to arrive at dynamical 
equations that can be obtained by employing contracted Bianchi identities. 
Conservation laws play signihcant part in establishment of instability range 
by more generic analytic approach, the dynamical equations in our case are 


■ 1 1 (^)' 

cyv^ = 0 ^ + —TO" (-R) = 0, 

Jr Jr 


1 {D) 

rrilv 


CyXu = 0 ^ {B) = 0, 

Jr Jr 


1 (o)' 

rTi2v 


G'‘:L„ = 0 ^ —R-'- + —R^ (rB) = 0, 
JR Jr 


(2.15) 


(2.16) 


(2.17) 


on simplihcation, we have 


GOo + goi + G0| + G° 


3A^ ^ 2B^ ^ G® 
^ ^ ^ "G 




2A 2B G \ 


3R' 2B' C 1 

“X ^ ^ "G ^ r 


^o.yBB _,oGG 

+ ^ ^ ^ = 0’ (2.18) 


A C AB 


A C B 


A' 3B' C 


^12 fA^ 4R® Cy ^22 f 

+r t+x+g 


ABC 


+ G00^ = 0, (2.19) 
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A' 

'a 


Go + G',! + G^2 + G 
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^2J^2 
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T ^ “b“ ^ "C 


B ' C 
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+ 


(^(^9 

- G ^^^2 = 0 - ( 2 - 20 ) 


The notations of 0,1 and 2 indicates t, r and 6. Terms belonging to matter or 
effective part of the dynamical equations can be viewed separately by insert¬ 
ing components of Einstein tensor given in Eqs. fl2.7p - fl2.13p . The dynamics 
of gravitating axial system can be explored with the help of perturbation 
scheme, that is useful in estimating the change in system with the passage 
of time. 


3 f{R, T) Model and Perturbation Approach 

Selection of the model under observation is a crucial constituent of the anal¬ 
ysis. Since we are dealing the system analytically, thus model selected 
shall bring fruitful mechanism for some particular form of f{R,T). We 
found that the f{R,T) form suitable for dynamical analysis is constrained 
to f{R,T) = f{R) -|- AT, where A is positive constant and f{R) is an arbi¬ 
trary function of Ricci scalar. The origin of such restriction comes from the 
fact that non-linear terms of trace in f{R,T) complicates the formation of 
modified field equations that can not be handled analytically. Such f{R,T) 
models bearing non-linear terms of trace of energy momentum can be dealt 
by using numerical techniques leading to more specified outcomes, whereas, 
the findings of analytic approach yield more generic. The viable f{R,T) 
model we have chosen is 


f{R,T) = R + aR‘^ + XT, (3.21) 

where any positive value can be assigned to a and A. 

The onset of modified field equations is non-linear in nature whose so¬ 
lution is still undetermined that is why perturbation approach is utilized to 
monitor variations in gravitating system with the time transition. All phys¬ 
ical quantities are taken to be time independent initially, passage of time 
implicates dependence on time as well. To introduce hrst order perturba- 
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tions, we chose 0 < e 1 


A(t,r,9) = Ao{r,9) + eD(t)a{r,9), ( 3 . 22 ) 

B{t,r,9) = Bo{r,9) + eD{t)b{r,9), (3.23) 

C(t,r,9) = Co(r,9) + eD(t)c(r,9), ( 3 . 24 ) 

p(t, r, 6) = pQ{r, 6) + ep{t, r, 6), ( 3 . 25 ) 

Pxx(j'')PxxO (r,9) + ePx^(t,r,9), ( 3 . 26 ) 

Pyy(t,r,9) = Pyyo(r,9) + ePyy(t,r,9), ( 3 . 27 ) 

Pzz(tA:^) = Pzzo(r,9) + eP^^(t,r,9), ( 3 . 28 ) 

Pxy(t,r,9) = Pxyo(r,9) + ePxy(t,r,9), ( 3 . 29 ) 

R(t,r,9) = RQ(r,9) + eD(t)e(r,9), ( 3 . 30 ) 

f{R,T) = [Ro{r,9) + aRl{r,9) + XTo{r,9)]+eD{t)e{r,9)[l 

+ 2aRo(r,9)], (3.31) 

fji = 1 + 2aRo(r,9) + e2aD(t)e(r,9), (3.32) 

/t = A. (3.33) 


The hrst order perturbed Bianchi identities fl2.18p - fl2.20p imply 
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(3.36) 


where Zip, Z 2 p and Z^p given in appendix. To make simplification a bit easier 
we substitute, / = 1 + 2aRo and J = e2aRo- The expression for energy- 
density p is derived from Eq. fl3.34p as 


P 



a 2Xih AicA 

-^- — + — 

Alo Bo Co) 


+ -J^{PxxO + Pyyo) + -p^PzzO 
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(3.37) 


The energy density and pressure stresses are associated as [28l [36] 

p. = r p. (3.38) 

PO + PiO 

where T describe the variation of different stresses with energy density. The 
index has variation as i = xx,yy,xy, zz, and Eq. fl3.37p together with Eq 
(I3.38P leads to corresponding perturbed stresses. Implementation of linear 
perturbation on Ricci scalar yields an ordinary differential equation having 
solution of following form 

D{t) = -e^\ (3.39) 

The expression for Z 4 is provided in appendix, Eq. fl3.39p is valid for overall 
positive values of Z 4 . 
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4 N and pN Approximation 

This section constitutes the terms belonging to N and pN limits with insta¬ 
bility criterion in the form of adiabatic index. Making use of Eqs. fl3.39p 
and (I3.38P in Eq. (13.351) leads to the evolution equation. The N and pN 
approximations for considered system are discussed in following subsections. 


4.1 Newtonian Approximation 

To approximate instability/stability range in Newtonian regime, we let Aq = 
1, Bq = 1, pq ^ piQ] i = XX, yy, xy, yy and Schwarzschild coordinates Cq = r, 
evolution equation along with these assumptions yield 


T < 


AiV' - fiVo - 2A(poiV2)' - l{P.yoN2y + AiVsiVg - fiV 2 + XP^yoN^ + Z^M 


XyPxxoNl)' + ^{PxyoNlY — \.Nl{PxxO + PyyO + Pzzo) 


where Z 2 N corresponds to the N-approximation terms of Z 2 p, and 


(4.40) 
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= a2Ai6, N2 = h + y^, 

r 1 


iVs = a' + 46' + 2(-)', N^ = a^ + Ab^ + - 


The inequality for T contains both material functions and effective part en¬ 
tries, system remains stable as long as the inequality fl4.40p holds. The terms 
appearing in expression for T are presumed in a way that all terms maintain 
positivity, this requirement impose some restrictions on physical parameters. 
The constraints in N-approximation are 


PxxO + PyyO + PzzO < -J^HPxxoNi)' + < ~2A(po-A^2)^ 

Violation of these constraints imply instability in the sources and thus lead 
to gravitational collapse. 
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Post Newtonian Approximation 

In pN approximation, we assume Aq = 1 — ^ and Bq = 1 + ^, corresponding 
inequality for range of stability is 


r < 


+ -^9-^10 ~ 2A(po-^6)^ ~ -{PxyoNo) + XpoN-j — -Nq + XP^-yoNg + Z^pN 
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(4.41) 
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Likewise Newtonian limit metric coefficients and effective part terms can be 
constrained to maintain stability of self gravitating system. System is stable 
unless above mentioned inequality holds, system collapses when ordering 
relation fl4.41l) breaksdown. One can deduce results of GR approximations 
by choosing vanishing values of A and a. 


5 Summary and Discussion 

Observational Signatures supports the argument that gravitating sources 
might deviate from spherical symmetry incidently. Thus non-spherical sym¬ 
metries facilitate in examining the realistic situations such as large scale 
structures, weak lensing, CMB etc. Motivating from the signihcance of non- 
spherical symmetries, we intend to explore impact of axially symmetric grav¬ 
itating source in context of f{R, T) gravity. More particularly we are dealing 
with restricted axial symmetry by ignoring meridional motions and rotation 
about symmetry axis. The consequence of restricted character of spacetime 
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leads to vorticity-free case, because absence of dtdO and dtdcj) terms indicates 
that vorticity of gravitating source vanishes for at rest observer. The met¬ 
ric under consideration is axially symmetric with three independent metric 
functions. 

Implications of axial symmetry on gravitating system has been studied ex¬ 
tensively in GR and modihed theories of gravity. The alternative gravity the¬ 
ory we have chosen to establish instability range is f{R,T) gravity, because 
dynamical instability of axially symmetric sources in f{R,T) framework has 
not been ascertained yet. The model under study f{R, T) = R + aR^ + XT 
is viable for positive values of a and A. The modihed held equations are 
obtained by varying action 03.281) for anisotropic matter distribution. The 
components of held equations 02.7I) - 02.13I) are used to arrive at conservation 
equations 02.18p - fl2.20p . These equations are of fundamental importance in 
establishment of instability range analytically. 

The held equations are non-linear in nature, its a difficult task to eval¬ 
uate their general solution. To count with this issue, we consider linear 
perturbation of usual matter and dark source terms. The perturbed phys¬ 
ical quantities such as energy density and anisotropic pressure stresses are 
extracted from linearly perturbed components of held equations, that are 
further inserted in perturbed Bianchi identities to arrive at collapse equation 
carrying both material and dark source ingredients. An ordinary diherential 
equation is formed from perturbed Ricci scalar, whose solution together with 
evolution equation provides adiabatic index. 

Adiabatic index dehnes range of instability for N and pN approximations 
inducing some constraints on physical quantities that are provided in previous 
section. Corrections to GR and f{R) gravity can be determined by setting 
a —)■ 0, A —)■ 0 and A —)■ 0 respectively. 
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Appendix 

Following equations contain linearly perturbed terms of conservation equa¬ 
tions and Ricci scalar respectively. 
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